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It is shown that one can in principle constrain the CP-violating parameter 5 from measurements 
of four independent |Vi,-| 2, s, or three \Vij\ 2 and a ratio of two of them, in the leptonic sector. To 
quantify our approach, using unitarity, we derive simple expressions in terms of four independent 
\Vij\ 2 's for cos 5, and an expression for sin 2 S from J 2 . Thus, depending on the values for \Vij\ and 
their accuracy, we can set meaningful limits on To illustrate numerically, if |V^,i| 2 is close to 
0.1 with a 10% precision, and if | V^3 1 2 is larger than 0.005 and for values of \V e 2\ 2 and | 1 2 that 
stay within ±0.1 of the current experimental data leads to a bound: n/2 < \8\ < it. Alternatively, a 
certain combination of parameters with values of | V^3 1 2 larger than 0.01 leads to a closed bound of 
73 < \5\ < 103. In general, we find that it is better to use |V^i| 2 or |K-i| 2 as the fourth independent 
\Vij\ 2 and that over most of the parameter space, S is least sensitive to |T4s| 2 . With just three 
independent measurements (solar, atmospheric and reactor), it is impossible to set limits on the 
CP phase. In this respect, we study the use of ultra high energy cosmic ray neutrino fluxes as the 
additional fourth information. We find that within the SM, neutrino fluxes of all three flavours will 
be very similar but that pushing current neutrino data to their extreme values still allowed, ratios of 
cosmic neutrino fluxes can differ by up to 20%; such large discrepancies could imply negligibly small 
CP-violation. We also study a non radiative neutrino decay model and find that the neutrino fluxes 
can differ by a factor of up to 3 within this model and that an accuracy of 10% on the neutrino 
fluxes is sufficient to set interesting limits on 5. 
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I. INTRODUCTION 

One of the significant achievements of particle physics 
has been our better understanding of the leptonic flavor 
mixings. Following some very detailed and pains taking 
measurements the neutrino anomaly has now been con- 
firmed; both in the solar and atmospheric sector [1]. An 
immeadiate conclusion is the existence of new physics be- 
yond standard model and to introduce a small neutrino 
mass. Such an extension leads to the notion of neutrino 
oscillation [2] similar in spirit to quark sector. A genuine 
three flavor analysis, allows for a solution space consist- 
ing of two mass-squared differences and three angles. The 
solution amounts to finding the allowed parameter space 
for the solar mixing angle, and the atmospheric mixing 
angle along with a strongly constrained reactor angle. 
Phenonemologically, best fit values and allowed range for 
these three mixing parameters dictates a pattern that is 
almost being negligible (reactor) to moderate (solar) to 
maximal (atmospheric). For a review on the analysis, we 
refer to [3]. 

Similar to the quark sector, one can have CP-violation 
in the leptonic sector due to massive neutrinos. Unlike 
the quark sector however, the CP violation in the lep- 
tonic sector can come due to both Dirac and Majorana 
phases in the mixing matrix, depending on the nature 
of the neutrinos [4] Clearly, assessing the Dirac or Majo- 
rana nature of the neutrinos will be an important goal in 
the future. In this paper, we consider only CP-violation 
a la Dirac and do not address at all Majorana phases. 



As a pertinent question, how to measure CP-violation in 
neutrino oscillations has attracted a lot of attention in 
the past [5, 6] and no doubt, the search for leptonic CP- 
violation will be one of the main goals of experimental 
particle physics in the years to come. Given the strong 
reactor constraints [7] it is very difficult for current neu- 
trino oscillation experiments to measure CP-violation. 
Nonetheless, as a proposal, a measurement is possible 
by one searching for differences between neutrino and 
anti-neutrino survival in a long-baseline experiment and 
measuring the spectrum [8] . 

In the present analysis, we want to explore the pos- 
sibility of extracting information on CP violation in the 
leptonic sector without performing a direct CP violation 
experiment. Clearly, such an experiment will have to be 
done eventually but in the near future, it might be eas- 
ier to measure individual lepton mixing matrix elements 
(Vij). Thus, it becomes interesting to see how informa- 
tion about these can be translated into information about 
CP violation in the leptonic sector and what level of pre- 
cision will be required on the |14, | 2 's to set interesting 
constraints on 5. 

In the next section, we will contrast two parametriza- 
tions that represent leptonic mixing; the usual 
parametrization in terms of mixing angles and phases [9] 
and another one constructed purely as moduli elements, 
where CP violation is not explicit but still present. Fol- 
lowing this, we will derive expressions that allow one to 
extract information on CP violation from measurements 
of four |Vy|'s or combinations of them. These relations 
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are used to calculate the precision on the measurements 
of the \ Vij | 2 's required in order to set constraints on S. We 
will then consider the use of UHCR neutrino fluxes to ob- 
tain the fourth |V^| 2 needed to set constraints on S [10]. 
We also derive general relations about their fluxes and 
estimate the precision required in order to set interesting 
limits on 6. To conclude, we will also consider a specific 
model that digresses from the SM and will strongly af- 
fect the neutrino fluxes that reach the earth and see how 
useful this can be in our extraction of the phase 5. 



II. LINKING FOUR \V tJ \ 2 AND 8 

A. Rephasing invariant parametrization 

We assume that we have three neutrino flavours and 
consider only CP-violation a la Dirac. We also assume 
unitarity in the mixing matrix even though some popular 
models of neutrino mass generation lead to violation of 
unitarity. 

The elements of the lepton mixing matrix are defined 

as 

v e \ fV el V e2 V e3 \ f Vl \ 
r>n ]=[ V V„ 2 V^ 3 \\V2 • (2-1) 

VrJ \V T 1 V T 2 V t3 ) \ V 3 I 

It has been shown before [11] that all the information 
about mixing and CP violation in the quark sector can 
be parametrized in terms of four independent moduli of 
(2.1). These moduli are rephasing invariant and basis in- 
dependent. A set of four independent parameters is not 
unique and any set of four parameters where we don't 
have three on the same row or column is acceptable. 
Thus, there are nine such ensembles which are allowed. 
For the present analysis, our choice is \V e2 \, [V^l, |V^i| 
and This choice seems appropriate since we have 

experimental information on three of the four parame- 
ters. 

The exact expressions between the moduli of the el- 
ements of the mixing matrix and the mixing angles are 
obtained to be: 

\V e3 \ 2 = sin 2 9 13 , (2.2) 
\V e2 \ 2 = sin 2 6 12 cos 2 6 13 , (2.3) 
|V^ 3 | 2 = sin 2 6 23 cos 2 9 13 . (2.4) 

Our current knowledge on the magnitude of the ele- 
ments of the leptonic mixing matrix comes from experi- 
ments on neutrino oscillation and can be summarized as 
follows [3]: 

sin 2 6»i2 = 0.3 ±0.08 
sin 2 6> 23 = 0.5 ±0.18 
sin 2 6»i3 < 0.05 at 3ct (2.5) 



B. Unitarity and cos(S) 

At the present time, there is no experimental infor- 
mation about |V^i and one can only limit its modulus 
through unitarity. In the notation of (2.1), we can write 

V,i = - 1 _ ^ Ve3l2 (\V e2 \\V T3 \ + \V el \\V, 3 \\V e3 \e i5 ) , 
\V e i\ = Vl-\Ve2\ 2 -\V e3 \ 2 , 

[Ksl = ^l-\V, 3 \ 2 -\V e3 \ 2 . (2.6) 
Expanding up to order | Ve3 1 2 , we calculate \V^\\ 2 to be 

|^!| 2 = |^ 2 | 2 (i-|^ 3 | 2 ) + |K 3 | 2 (|v;2| 2 -3|v; 2 | 2 |y Al 3| 2 

+ |T^3| 2 ) + 2|V e3 ||T42||^3| 

x cos^l-l^slVl-l^l 2 (2.7) 

With current experimental data (2.5) and adding the er- 
rors in quadrature, the limits on |V^i| 2 translate to 

IV^iI 2 = 0.15 ± 0.07 + (0.35 ± 0.15)|K 3 | 2 

± (0.46 ±0.03) |K 3 1 . (2.8) 

If we use the 3 a limit of 0.05 on \V e3 1 2 , we obtain a lower 
limit of and an upper limit of 0.35 on |V^i| 2 while the 
1 a limit of 0.012 on |V e3 | 2 leads to 0.029 < \ 2 < 0.28. 
Therefore, we can arrive at closed bound 

0.17 < \V^\ < 0.52 at la . (2.9) 

We can use this expression to set a bound on | V M i | 2 or we 
can invert it to get information on cos(S) once we have 
information on all four \Vij\ 2 , the exact relation being: 

, l^il 2 HK 2 | 2 (i-|^3l 2 ) + ff \v e3 \ 2 + c |y e3 | 4 

coso = — — — , 

x = 2\V e2 \\V e3 \\V^ 3 Wl-\V e2 \ 2 -\V e3 \ 2 

x ^i-\v e3 \ 2 -\v, 3 \ 2 , 
b = \v e2 \ 2 -\v^ 3 \ 2 -2\v^\ 2 + \v e2 \ 2 \v^\ 2 , 

c = |V,i| 2 + |^ 3 | 2 . 

(2.10) 

Clearly, without information on |V^i| 2 , we cannot say 
anything on 5. 

C. CP- Violation and the |V^|'s 

The measure of CP violation is expressed, in general, 
through the Jarlskog parameter defined as [12] 

J = %m(V e2 V„ 3 V; 3 V; 2 ) (2.11) 

Once we have chosen our four independent parame- 
ters of the mixing matrix, J 2 is not independent but an 
explicit function of these four parameters, namely [11] 
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have, we are infact setting sin 2 (S) = 0.25. The straight 
line in the middle is the curve \V ll i\ 2 j2=J 2 and is given 
(1- l^esl 2 - |V e2 H|K3| 2 |^i| 2 |Vv3| 2 by 2.7 with cos (S) = 0. 



- - 4 (\v e2 \ 2 \V^\ 2 + |v^| 2 |v; 3 | 2 - |u, 3 | 2 |u e3 | 2 



2, T/ 2^2 rolo ^ in. CONSTRAINING 5 FROM THE | 2 OR 



K 3 \'\V e2 \ 2 y . (2.12) 



THEIR RATIOS 



This is a general expression based on rephasing invari- 
ants and we can use it to set limits on |V M i| 2 : requiring 
positivity of J 2 leads to (2.7) with cos 5 replaced by ±1. 
We can also calculate what value of |V^i| 2 will maxi- 
mize J 2 such that |V^i| 2 2= j2 • Differentiating J 2 with 
respect to |V^i| 2 and retaining only first order in |V e3 | 2 , 
one then recovers the first two terms of (2.7). This is con- 
sistent since maximizing J 2 requires sin 2 5=1. One ob- 
serves here that a value of | V^i | 2 larger than | V^i | 2 2=J 2 
requires <5 to be in the first or in the fourth quadrant. 
This simple point does not require particularly precise 
data and could already be useful information for model 
builders. 

In the standard basis, J 2 is given by [13] 

J 2 = sin 2 {0 l2 )cos 2 {9 l2 ) sin 2 {d 23 )cos 2 (8 23 )sin 2 (9 13 ) 
cos 2 {9 13 ) cos 2 {9 13 )sin 2 {S) (2.13) 

Clearly, whether in the Vpmns or in our current basis, 
J 2 is the same physical observable and we can compare 
the two expressions for this parameter. Using the re- 
lations between the Vij and sin(8ij), we express J 2 in 
terms of sin 2 (6) and three of our |Vy| 2 : 

J2 = (1 _ |y e3 |2)2 l^ 2 | 2 |^ 3 | 2 |V/ e3 | 2 (l - \Ve 3 \ 2 - \V e2 \ 2 ) 

x (l-\V e3 \ 2 -\V^ 3 \ 2 ) S in 2 5 . (2.14) 

Using (2.12) and (2.14) we now have a relation between 
sin 2 S and our four \Vij\ 2 . Therefore, we can relate S to 
our set of four \Vij\ 2 through either cos S or sin 2 8. Since 
we want ultimately information on S itself and there is 
more information through cos 6 than through sin 2 (S), in 
what follows, we will concentrate on cos S. We will be 
left with a twofold degeneracy as we will not be able to 
get the sign of S. 

On figure 1, we plot J 2 as a function of \V^i\ 2 and 
|V e3 | 2 for central values of | V^ 3 1 2 (0.5) and \V e2 \ 2 (0.3). 
We see clearly that the vanishing of |U e3 | 2 leads to the 
vanishing of J 2 , as it must be. We also see that J 2 is a 
fairly linear function of |U e3 | 2 ; this is expected since the 
maximum value |U e3 | 2 can have is rather small compared 
to the other parameters and an expansion to first order 
in |U e3 | 2 would be adequate. For numerical purposes, we 
note that the largest value taken by J 2 on this figure is 
2.42 x 10~ 3 ; it is clearly at the largest possible value of 
|U e3 | 2 and at \V^\ 2 -0.167 

On figure 2, we plot curves of constant sin 2 (5). When 
we fix Vg2 1 2 , |U m3 | 2 , and V^3 1 2 , if we ask what value of 
\V^i\ 2 will give 25% of the maximum value that J 2 can 



In order to be able to set some bounds on 5, we need 
at least lower limits on |U e3 | 2 and |V^i| 2 besides the data 
that we have on \ V e2 \ 2 and |V^ 3 | 2 ; as we can see on figure 
2, just an upper limit cannot constrain S. In what follows, 
we explore the full 4-dimensional parameter space: we al- 
low \V^i\ 2 to cover the whole range allowed by unitarity, 
and |U e3 | 2 varies from 0.001 to 0.03, since 0.03 is close to 
the current 3-cx limit and 0.001 makes CP-violation ex- 
tremely small at planned CP-violating experiments[14]. 
As for \V e2 \ 2 and |V m3 | 2 , we will work with values that 
can differ by ±0.1 from their current central values of 0.3 
and 0.5, respectively. 

We assume that we have four |Vy| 2 with their ex- 
perimental uncertainties. The experimental central val- 
ues of the |Vy| 2 's (|Vij| 2 ) lead to the central value 
of cos S. In order to estimate the error, or range 
that we should associate to this central value, we 
use a Monte Carlo technique and we cover, for each 
\Vij\ 2 , the space (|Vy| 2 — experimental error) - (|Vy| 2 + 
experimental error). Within this 4-dimensional space, we 
are interested only in the largest and smallest values that 
cos S can have; these become the range that cos 5 covers 
with this particular combination of |Vy | 2 and their asso- 
ciated errors. 

As our fourth \Vij\ 2 , we picked | | 2 . This choice is 
not unique and we could have chosen |V^ 2 | 2 , or |Uri| 2 , 
° r \Vt 2 \ 2 - It is straightforward to rewrite all our equa- 
tions in terms of these parameters through the following 
relations: 

|U, 2 | 2 = l-|U Ml | 2 -|U, 3 | 2 

|Ki| 2 - \V e2 \ 2 + \V e3 \ 2 -\V^\ 2 

\V T2 \ 2 = \V^ 3 \ 2 + \V^\ 2 -\V e2 \ 2 (3.15) 

We studied all four parameters and we can say that, 
when they (ie |U M i| 2 , |U m2 | 2 , |Ki| 2 , IK2I 2 ) have the same 
experimental error (eg 5%) then: 

• if the known parameters stay close to their current 
central values, in general, |V^i| 2 leads to slightly better 
bounds than \V T i\ 2 

• if the known parameters digress substantially from 
their current experimental values, then \V T \\ 2 leads to 
slightly better bounds than iV^il 2 

• I V^j,2 1 2 and |V^- 2 | 2 are not as good as the previous two: 
the limits obtained from these parameters are degraded 
by about a factor of 2 when compared to those obtained 
from the previous two parameters; |V Al2 1 2 is, in general, 
a little bit better than |Ur 2 | 2 

• things improve a bit for the last two parameters if we 
give all four parameters the same absolute uncertainty 
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instead of the same relative uncertainty (eg we compare 
|V^i| 2 = 0.150 ±0.005 and \V t2 \ 2 = 0.350 ±0.005 instead 
of ±5% for all parameters). Even then, |V^i| 2 and \V T i\ 2 
remain better than the other two parameters, but by a 
factor of 1.5 instead of a factor of 2. 

We also take into consideration some ratios of 
\Vij\ 2 as potential fourth parameter. The first one is 
Pi = \V T i\ 2 /\Vfj,i\ 2 and has very similar properties to 
|K 2 |7I^2| 2 ; the second one is p 2 = \V t2 \ 2 /\V Tl \ 2 and 
has very similar properties to | 1 2 / 1 V^x | 2 - We also 
studied |U t2 | 2 /|U m i| 2 but it turned out to be too sensitive 
to both | V^=3 1 2 and | 1 2 to be of any use; and similarly 
for \V T i\ 2 /\V^ 2 \ 2 . The same can be said of any ratio 
that involves |V e i| 2 . So, we present in Tables I, and II 
the limits that we can set on cos<5 from measurements 
on U M i| 2 , p\ and p 2 . These are representative of what 
can be achieved. After studying the parameter space 
described above, we can say that: 

• uncertainties of 10% on the four parameters can lead 
to very tight bounds on S; the uncertainty on | Ve3 1 2 can 
be much larger than this without affecting the bounds 
very much. 

• an interesting constraint already occurs for the com- 
bination (0.3,0.5,P,|V e 3| 2 ) with a 10% uncertainty on 
the parameters: the first two parameters (|V e 2| 2 and 
| 1 2 ) are at their current central values while the third 
one (either |U M i| 2 = 0.15, p\ = 1.0 or p 2 = 2.3) is close 
to the value that maximizes J 2 . We find that if | V^3 1 2 
turns out to be large (0.03), then \S\ has a range of 40-70 
degrees centered at about 90 degrees, depending what 
P one uses; this range decreases as 1 1/^3 1 2 increases. If 
the uncertainties are reduced to 5%, the range becomes 
30-60 degrees and | V^3 1 2 can be reduced to 0.01. 

• for most combinations of \V e2 \ 2 and | V^3 1 2 , if |V^i| 2 
turns out to be relatively small, about 0.1, then \S\ has 
to be between ir/2 and it for any | V^3 1 2 larger than 
0.005. 

• the limits are not very sensitive to the uncertainty on 
| V^3 1 2 ; going from 5% to 20% does not change the limits 
by much. If the other uncertainties are small (2%), the 
uncertainty on | V^3 1 2 can go up to 50% and it is still 
possible to get interesting bounds on 5 

• the limits are not very sensitive to |U e 3| 2 itself and 
going from 0.01 to 0.005 will not change the limits very 
much. When \V e 3\ 2 becomes 0.001, however, the limits 
are somewhat degraded but still useful. 

• as the value of | V^3 1 2 becomes smaller, the uncertainties 
on the other parameteres must decrease in order to keep 
interesting bounds on 5; note that the uncertainty on 
| V^3 1 2 can be rather large without affecting much the 
bounds 

• in general, we find that the best parameter is p\, the 
second best is p 2 and the third best is |V^i| 2 . 

The previous analysis is very general and describes 
which parameters must be measured and with what 
accuracy in order to set limits on S. We now turn to 
some potential processes that could give us the fourth 



information that we need in order to be able to set limits 
on 5. 



IV. ULTRA HIGH ENERGY COSMIC RAYS 

Ultra high energy cosmic rays and their detectors 
have attracted a fair amount of attention in recent 
years. [15, 16, 17, 18] In particular, the use of UHCRs 
with the aim of observing leptonic CP-violation has been 
considered in [? ] An interesting consequence for the 
UHCR neutrino spectra is due to maximal atmospheric 
mixings. Following maximal mixing, or tan 6 a ~ 1 leads 
to a unique prediction for UHCR neutrino fluxes. It 
was shown that UHCR neutrinos (which are expected 
to be sourced by cosmic objects such as AGNs) when 
measured by ground based detectors, the expected fla- 
vor ratio <j> e : cf)^ : <j> T = 1 : 1 : 1 [20]. This value 
is also known as the standard flavor ratio and by itself 
constitute an independent confirmation of the neutrino 
mixing data from atmopheric sources. Since we propose 
to extract the Dirac CP phase present in the conven- 
tional PMNS mixing matrix, our discussion will be re- 
stricted to Dirac neutrinos. Before proceeding with the 
analysis, let us briefly allude to the significance of this 
proposal. It is well known that CP violation (due to os- 
cillation) in the leptonic mixings (even if larger than the 
quark secotr) will nonetheless remain a hard problem to 
resolve[14]. This arises from the strong constraints which 
reactor neutrinos set on the mixing element V e z [7]. As 
a consequence of this, even if the CP phase is large, the 
smallness of V e z leads to phase being insensitive to any 
CP violating measurements. 

A. The basic formalism 

In the same fashion as quarks mix, massive neutrinos 
mix (and also oscillate) between two eigenbas2s, the fla- 
vor (v a ) and mass eigenbases (vi). The two bases are 
related by a unitary matrix V such that v a = V a iVi- 
Here we assume the summation over the mass eigen- 
states. Corresponding to a particular mass eigensatcs 
is a mass value, rrii which in the limit of small mixings 
determines the mixing angles. Thus, in the limit of small 
mixings, 6s mixes eigenstates v\ and v 2 , while 6a mixes 
eigenstates v 2 and ^3 and 6r mixes the eigenstates V\ 
and z/3. 

In the context of UHCR neutrinos which travel as- 
tronomical distances, the coherence between the various 
mass eigenstates is averaged out. In other words, after 
production these neutrinos essentially travel as individ- 
ual mass eigenstate. At the point of detection, they are 
in the flavor states. Therefore, the detection probability 
in a given flavor eigenstate is 

4> e = 1 + 2x(2c 2 A - 1) ; x = (s s c s ) 2 , 
</y = 2xc 2 A + 2(c A (l-2x) + s\) , 
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4> T = s\ + 2x^(1 - c\) , (4.16) 

where s and c denote sine and cosine, respectively. From 
the above expressions and the experimental fact 8a — 
7r/4, it follows that all neutrino flavors must be detected 
with the same weight factor. In deriving this result, we 
have disregarded the mixing corresponding to reactor ex- 
periments, which is consistent with zero [7]. This result 
is also independent of the solar angle 6s- 

B. Neutrino Fluxes and \Vij\ 2 

Let us revise in more detail the arguments from the 
previous section. Consider the probability of oscillating 
from flavor a to (3 

3 

i>j=i 

3 

+ 4 ^ lm\K a(i ^\ sin(Aij) cos(Aij) , 
i>j=i 

K aPtij = VaiV&V^Vn ; A y = % 3 . (4.17) 

We can express the real and imaginary parts of K a p t ij 
in terms of moduli as follows: 

2Re[K a/3 ^] = \V ai \ 2 \V 0j \ 2 + \V fii \ 2 \V aj \ 2 

- ^2c af}7 c ijk \V 7k \ 2 (4.18) 

and 

Im[K a 0^j} = J e a p Sij (4.19) 

where we have defined: 

Cijk = 1 if i j, j k, k^i 

= otherwise (4.20) 

and e as the antisymmetric tensor 

/ -1 1\ 
e = 1 0-1 (4.21) 
V-l 1 0/ 

We are interested here only in the probabilities for 
v t-, Vfj, — > v e and v T — ► v e . After combining two 
terms (A13 = A23 to high accuracy), expressing the third 
term as a combination of |Vy| 2 and averaging over the 
oscillating terms, we obtain 

Put = 2\V^ 3 \ 2 \V T3 \ 2 ~ (\V el \ 2 \V e2 \ 2 ~ IV^IV^I 2 

\Vn\ 2 \V T2 \ 2 ) , 
P eT = 2|t/ e3 | 2 |T/ r3 | 2 -(|V- All | 2 |y Al2 | 2 -|T4 1 | 2 |K 2 | 2 

- \VtA 2 \Vt2\ 2 ) , 

Pe» = 2|^ 3 | 2 |K3| 2 -(|Kl| 2 |y r2 | 2 -|V| 2 |^2| 2 

- \V el \ 2 \V e2 \ 2 ) . (4.22) 



One then uses the relations 

\v el \ 2 = i-\v e2 \ 2 -\v e3 \ 2 , 
\v, 2 \ 2 = i-Iv^i 2 -!^! 2 , 
IK1I 2 - \v e2 \ 2 + \v e3 \ 2 -\v^\ 2 , 
IK2I 2 - \v lil \ 2 + \v ll3 \ 2 -\v e2 \ 2 . 

(4.23) 

to express these probabilities in terms of our four |Vy | 2 . 

Note that when averaging over the oscillating terms, 
the third term of 4.17 vanishes; therefore, one will not be 
able to observe directly CP-violation with UHCR's. 

If we assume that the main production mode of these 
UHCR is [21] 

p + 7-»7r ± +X , (4.24) 

with subsequent decay of the n ± into muons, electrons, 
and neutrinos, we conclude that there will be two (or 
9^) for every v e (or v e ) and virtually no v T (or v T ). The 
initial fluxes are then cfP e , <f>° = 2<fi®, and <fi® = 0. In 
order to calculate a given neutrino flux that reaches the 
earth, we take into account the probability that this neu- 
trino will oscillate into other types of neutrinos and the 
probabilities that other neutrinos will oscillate into this 
type of neutrino. Up to a common geometrical factor, 
the observed terrestrial fluxes are: 

4>* e = 4>l (l+Pe^-Per) , 

4>l = <£ (2 - P e , - 2P, T ) , 

# - 4>l (2P^ + Per) . (4.25) 

In terms of matrix elements, we get 

& = 4P e (l-2\V e2 \ 2 {\V, 3 \ 2 -\V e2 \ 2 ) 

+ 2|Vv 1 | 2 (l-2|V/ e2 | 2 ) 

- 2|V/ e3 | 2 (l + |V^i| 2 - IK2I 2 - \V^\ 2 ) + 2|V e3 | 4 ) , 
<j>l = 4P e (2 + (l-|Vv 3 | 2 )(|V/ e2 | 2 -4|V, 3 | 2 ) 

- |^i| 2 (3 + 2|F e2 | 2 -4|Vv 3 | 2 -4|Vv 1 | 2 ) 

+ |K3| 2 (|^| 2 -|VI 2 )) . 

£ = <£ (4|^ 3 | 2 (i-|V M3 | 2 )-|K 2 | 2 (i-|K 2 | 2 ) 

+ 3|y e2 | 2 (|^ 3 | 2 )-|K 2 | 2 ) 

+ |K 3 | 2 (2 + 3|^!| 2 - 3|^ 3 | 2 - 2|Fe 2 | 2 - 2\V e3 \ 2 ) 

+ |U All | 2 (l-4|Vv 1 | 2 -4|y M3 | 2 + 6|V; 2 | 2 )) (4.26) 

Clearly, we don't know the original neutrino fluxes at 
the source. Therefore, it is more meaningful to try to 
measure ratios of neutrino fluxes. We will consider the 
ratios R efl = </>'/<^i an d Per — 4^1 4> r an d calculate how 
they vary when we cover the parameter space available. 
Our results are summarized in Table III. One can see 
that: 

• by numerical accident, |V^ 3 | 2 close to 0.5 has a very 
strong influence on R efl and R eT and tends to keep them 
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close to 1 

• if | V^3 1 2 stays close 0.5 to within 10%, R efx and R eT 
remain close to 1, for any value of \V e 2\ 2 and \V 2 3 \ that 
are within the experimentally allowed ranges 

• when one varies |V^3| 2 by 20%, either up or down, it 
is possible for R efl and R eT to differ from 1 by as much 
as 20%. Unfortunately, these large discrepancies require 
values of iV^il 2 such that S ~ or n, which means very 
little CP-violation 

• | V e3 | 2 has very little impact and varying it from 0.01 
to 0.005 will change R efl and R eT , typically by 1% or so. 



Therefore, we can say that even within the stan- 
dard production mechanism for the UHCR, deviations 
of R eix and R eT by 20% from their expectation value of 
1 are possible with the current experimental data but 
require in general pushing these data to their limits and 
CP-violation to be vanishingly small. 



C. Neutrino Fluxes and S 

In order to study the bounds one could set on 6 if one 
were to measure R e ^ with a given accuracy, we express 
|V^i| 2 in terms of R efl and use that result directly into 
our expression of cos(5) (2.10). Again, we use a Monte 
Carlo technique to scan the allowed parameter space. 
As a technical detail, we note that we have a quadratic 
equation to solve with two possible values for |V^i| 2 
but only one value respects unitarity and leads to a 
positive value of J 2 ; events that do not respect unitarity 
are excluded from our Monte Carlo. In Table IV, we 
summarise our results. We see that: 

• the most important parameters are |V^3| 2 and R efl , 
where the uncertainties must be rather small, while 
| V^2 1 2 and | V^3 1 2 can tolerate much larger errors without 
degrading the limits much 

• in general, a precision of two percent on all 4 parame- 
ters leads to tight constraints on 6; this can be relaxed to 
10% for \V e 2\ 2 and up to 25% for \V e3 \ 2 without loosing 
much on the bounds on 6 

• large or small values of R eil (1.13 or 0.84, for example) 
require higher precision in the four parameters in order 
to keep cos(5) away from +1 or -1; the precision required 
on | Ve2 1 2 and | V^3 1 2 is still lower than that on the other 
two parameters. Recall that cos(S) = ±1 means no 
CP-violation. 

• as | V^3 1 2 gets smaller, a higher precision on the 
parameters is necessary in order to keep the same tight 
range on 6 

• with some combinations of parameters, it is possible 
to exclude S = n or with a precision of a few percent 
on the parameters 

• an uncertainty of 5% on R efl can limit 5 to a 90-degree 
range or smaller with certain combination of parameters 



The precision required in order to get interesting 
constraints on S is daunting but this was expected since 
we saw that |V^i| 2 was the least sensitive of our three 
parameters to set limits on 5. 



V. NEW PHYSICS AND 6 

This general behaviour of cosmic neutrinos fluxes has 
been noted before [20] and comes directly from the fact 
that | V^j,3 1 2 is close to 0.5. As | V^3 1 2 goes away from 0.5, 
(current data allow a deviation of up to 20%) the ra- 
tios can differ substantially from 1:1:1, as we just saw. 
This brings the question of how one would interpret ra- 
tios that would differ much more than 8%, say 30% or 
more. It would be difficult to explain these anomalies 
within the SM and such discrepancies might suggest some 
new physics. This new physics could bring its own CP- 
violation, which could be a la Dime; we would then deal 
with an effective CP-phase such that S e ff ec ti ve = S+S new . 

We will consider new physics beyond massive neutri- 
nos and for illustration purposes, we choose non-radiative 
neutrino decays. 



A. Non-radiative neutrino decays and Seff 

The notion of non-radiative neutrino decay as a source 
of deviations from the standard ratio was first discussed 
in [22] and we refer to this paper for details. Essentially, 
in this model, the observable flavor ratio is modified to 
be 



= \V el \ 2 :\V, 1 \ 2 :\V rl \ 



(5.27) 



where we assume normal mass hierarchy. 

Contrary to the SM flux ratios, the ratio R T ^ = <j> T /<t>iJ. 
is a direct measure of 



Pi = 



\Vn\ 
l^i I 



(5.28) 



which proved to be the best variable to set constraints on 
5. Table V shows what kind of limits one could obtain 
with such a scenario. We note that: 

• uncertainties in the 5% - 10% range can lead to inter- 
esting constraints on <5 e //. 

• the limits are insensitive to the uncertainty on | X^;3 1 2 

• in this model, flux ratios of up to 3 are allowed with 
current data. This limit depends on | V^3 1 2 and as | V^3 1 2 
decreases, so does the upper limit of the ratio. 

As we just saw, if R e ^ or R er is close to 1, there is 
no need for new physics and a precision of 1% - 2% is 
required to set bounds on 8. Such a precision appears 
extremely difficult at planned detectors such as IceCubc, 
Antares, Nestor, Anita, or Baikal [23]. An optimistic 
view would be the detection of a few tens of events per 
year at IceCube, for example. If we assume that all these 
detectors will have comparable performance after a few 
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years of running, then we could hope for about 100 events 
per year. This will not be sufficient to reach a precision 
of 1% in a reasonable time frame but a 10% precision ap- 
pears within reach. Therefore, the better scenario, when 
considering the experimental precision, is to have very 
different neutrino flux ratios; in that case, a 10% preci- 
sion leads to interesting constraints on S e f fective ■ 

VI. CONCLUSIONS 

In this paper we have shown how it is possible to 
get information about CP-violation in the leptonic sec- 
tor without performing any direct experiment on CP- 
violation. In the usual parametrisation of the leptonic 
mixing, the mixings are described by three mixing angles 
and the fourth parameter is the CP-violating phase, S. 
All the information about mixing and CP-violation can 
also be described completely by four mixing elements, 
four | Vij | 2 . We have shown how one can extract infor- 
mation on sin 2 {5) and cos(S) from four \Vij\ 2 or three 
\Vij\ 2 and a ratio of two \Vij\ 2 . We used two such ratios: 
Pi = \V Tl \ 2 /\V^\ 2 and p 2 = \V T2 \ 2 /\V Tl \ 2 . 

We find that if one uses | V^ 2 1 2 , V^3 1 2 and |V e3 | 2 , then 
it is better to use |V^i| 2 or \V T \\ 2 as the fourth inde- 
pendent |Vy| 2 . When the first two parameters remain 
close to their current central experimental values, |V^i| 2 
is slightly better than \V T \\ 2 while \V T \\ 2 is slightly better 
when the first two parameters digress from their current 
central experimental values. 

As a first constraint that does not require particu- 
larly precise data, we find that if |V^i| 2 were larger than 
| V^ti 1 2 2 =J 2 , then S would have to be in the first or 
fourth quadrant. 

In general, when the four parameters have uncertain- 
ties of 10%, it is possible to get interesting constraints 
on S on some of the parameter space. For example, as- 
suming 10% errors, if | Ve2 1 2 and V^] 2 keep their current 
experimental central values and \V^i\ 2 turns out to be 
about 0.15, then \S\ has a 70 degree range centered at 
about 90 degrees if | V^3 1 2 turns out to be 0.03; this range 
decreases to about 40 degrees if one measures p\ and it 
decreases as |Vg| increases. If the uncertainties are re- 
duced to 5%, then the ranges on S become 30-60 degrees 



and can tolerate a value of 0.01 for IV^I 2 . We also found 
that, for most combinations of \V e 2\ 2 and |V^3| 2 , a small 
value of IV^il 2 (0.1) implies that \S\ is between n/2 and 
7r. It turns out that cos(S) is not particularly sensitive 
to | V^3 1 2 over most of the parameter space. This implies 
that once the errors are reduced to the 2% level, even if 
the error on |V e 3| 2 stayed at 10%, it would not degrade 
the constraints on 5 by much. Of the three variables 
investigated here, the best one turned out to be p\. 

In order to get the fourth piece of information needed 
to set bounds on 5, we explored the use of ultra high 
energy cosmic neutrinos. Generally, due a numerical ac- 
cident stemming mostly from |V^3| 2 ~ 0.5, the cosmic 
neutrino fluxes are very similar once they reach the earth. 
However, current data allow cosmic neutrino fluxes to dif- 
fer from their central, expected value of 1 by up to 20%: 
this requires pushing the parameters to the acceptable 
limits and implies that CP-violation in the leptonic sec- 
tor is vanishingly small. Unfortunately, the uncertainties 
on the data would have to be at the few percent level on 
the cosmic neutrino fluxes and on | V e 2 1 2 in order to set in- 
teresting limits on S. Finally, we studied a non-radiative 
decay model of neutrinos that could lead to flux ratios 
very different from 1. Indeed, in this model, the current 
data allow for neutrino fluxes to differ by a factor of up 
to three. Fortunately, in this model, one ratio of neutrino 
fluxes is p\ and leads to rather interesting constraints on 
5 e ff. with experimental uncertainties of 10% on ratios of 
cosmic neutrino fluxes. 

It appears that 1% - 2% precision on extremely high 
energy neutrino fluxes will be out of reach for the planned 
detectors but that 10% precision will be within reach. 
It would then be difficult to constrain S within the SM 
context using these cosmic rays and on would have to 
rely on other means to measure |V^i| 2 or | Wi | 2 or pi or 
p 2 - On the other hand, if the neutrino fluxes are very 
different, then one is likely outside the SM and a 10% 
precision on the neutrino flux ratios is sufficient in some 
models to set interesting bounds on 5 e f f ec tive ■ 
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Figure Captions 

Figure 1 

J 2 as a function of \V e 3\ 2 and |V^i| 2 for the current cen- 
tral values of \V e2 \ 2 (0.3) and |T^ 3 | 2 (0.5) 
Figure 2 

Different values of sin 2 {8) in the |V e 3| 2 - | V^i | 2 



plane. ^From the outside, inward, we plot sin 2 {8) — 
0,0.1,0.25,0.5,0.75,0.95,1. The line in the middle cor- 
responds to sin 2 (8) — 1 and is given by eq. 2.7 with 
cos(8) = 0; the most outward curves correspond to 
cos(8) = ±1. 



Table I 



\v e 2? 


V,s\ 2 


IK3I 2 


Pi, Pi, \V^i\ 2 


cos(5)- 


cos(5) c 


cos(5) + 


0.30 (10) 


0.50 (10) 


0.03 (10) 


pi 


= 1.0 (10) 


-0.257 


0.056 


0.387 








P2 


= 2.3 (10) 


-0.300 


0.111 


0.457 












2 = 0.15 (10) 


-0.708 


-0.132 


0.442 


0.30 (5) 


0.50 (5) 


0.01 (5) 


Pi 


= 1.0 (5) 


-0.224 


0.033 


0.294 








P2 


= 2.3 (5) 


-0.285 


0.042 


0.329 












2 = 0.15 (5) 


-0.569 


-0.076 


0.414 


0.30 (5) 


0.50 (5) 


0.01 (20) 


Pi 


= 1.0 (5) 


-0.247 


0.033 


0.310 








P2 


= 2.3 (5) 


-0.303 


0.042 


0.343 












2 = 0.15 (5) 


-0.607 


-0.076 


0.457 


0.30 (5) 


0.50 (5) 


0.01 (5) 


Pi 


= 1.0 (10) 


-0.306 


0.033 


0.391 








P2 


= 2.3 (10) 


-0.404 


0.042 


0.440 












2 = 0.15 (10) 


-0.747 


-0.076 


0.599 


0.30 (5) 


0.50 (5) 


0.015 (5) 


Pi 


= 1.0 (5) 


-0.170 


0.040 


0.255 








P2 


= 2.3 (5) 


-0.190 


0.065 


0.300 












2 = 0.15 (5) 


-0.493 


-0.093 


0.306 


0.30 (5) 


0.50 (5) 


0.005 (5) 


Pi 


= 1.0 (5) 


-0.339 


0.023 


0.392 








P2 


= 2.3 (5) 


-0.429 


0.007 


0.412 












2 = 0.15 (5) 


-0.756 


-0.054 


0.642 


0.30 (5) 


0.50 (2) 


0.001 (25) 


Pi 


= 1.0 (5) 


-0.529 


0.010 


0.555 








P2 


= 2.3 (5) 


-1.0 


-0.081 


0.885 












2 = 0.15 (5) 


-1.0 


-0.024 


1.0 


0.30 (2) 


0.50 (5) 


0.001 (25) 


Pi 


= 1.0 (5) 


-0.879 


0.010 


0.901 








P2 


= 2.3 (5) 


-0.796 


-0.081 


0.552 










V»x 


2 = 0.15 (5) 


-1.0 


-0.024 


1.0 



Values of cos(S) for different combinations of parameters: cos(5) c is obtained from the 
central values of the parameters while cos(5)- and cos(5)+ are the extremum values 
allowed by the uncertainties (shown in percent, in parentheses next to each parameter). 



Table II 



IK2I 2 

1 c:^ 1 


IK3I 2 

1 1 


IK3I 2 

1 1 


jOi , £>9, 1 V„i 2 


cos(5)— 


cos(5)r- 


cos(t5)_|_ 


0.30 (5) 


0.40 (5) 


0.01 (5) 


Pl = 0.68 (5) 


-0.204 


0.025 


0.247 








rtn — 9 9 


-n 9Q4 


-n n^4 


n iQ5 








\Y , I 2 — n 1« 

I ^/Ul| — VJ.-LO ^tJJ 


-n fii 5 


-n n7fi 


n 458 


n ^n (~\ n^i 

U.OU I J-VJ 1 




vj.vj ± \o j 


„ _ f) CO fKA 
pi — VJ.VJO I O J 


n 91 ^ 

vj.z ±0 


n n9^ 


n 9^4 

VJ.ZtJ^i 








^2 — ^-^ v,<-V 


-n 448 


-0 n^4 


n "^^n 

VJ.OOVJ 








IT/ |2 _ n 10 

| '/xl | — VJ.-LO ) 


-n 8i 3 


-n n7fi 

VJ . W I 


n 665 


n ^n 


VJ.4:VJ I J.VJ 1 


vj.vj-l \ 01 


px — VJ.UO V (J I 


n ^4Q 


n n9 K i 

VJ.VJZU 


VJ.OOO 








rtn — 9 9 l^l 

P2 — ^-^ v,<JJ 


-n ^1 n 

vj.o iu 


-n n^4 


n 1QQ 








IT/ 1 |2 _ n 10 /c\ 
1 v fil\ — vj.±o yo ) 


-n 7fi5 


-n n7fi 


n 582 


n ^n 


VJ.^iVJ Id J 


VJ. VJ _L I 1. VJ 1 


fj\ — VJ.UO 1 1 


919 


n n9^ 


n 9^n 

VJ.ZtJVJ 








rto — 9 9 


-n ^nfi 

VJ.OVJU 


-0 n^4 


n iQfi 

VJ. J. C/U 








\y 12 _ r. 10 (c\ 

1 ^/xl | — VJ.-LO ^UJ 


-n 628 


-n n7fi 


n 471 

VJ.t: f -L 


D 9^ (9~\ 


n ^ (9^ 

VJ. UtJ I Z 1 


n ni (9^ 

VJ.VJ-L \^ ) 


rt , - (IS (0\ 
Pi — u.o i^zj 


\J.UO J. 


U.UtJW 


VJ. 1 UVJ 








= 3 85 C2 1 


0.541 


0.652 


0.770 








|F M i| 2 = 0.145 (2) 


0.481 


0.663 


0.856 


0.25 (2) 


0.55 (2) 


0.01 (25) 


Pi = 0.8 (2) 


0.553 


0.650 


0.854 








p 2 = 3.85 (2) 


0.543 


0.652 


0.852 








|F Ml | 2 = 0.145 (2) 


0.482 


0.663 


0.993 



Values of cos(5) for different combinations of parameters: cos(S) c is obtained from the 
central values of the parameters while cos(5)- and cos(5)+ are the extremum values 
allowed by the uncertainties (shown in percent, in parentheses next to each parameter). 



Table III 



O 

o 

<D 

in 



(N 

> 

vo 
o 
in 
o 

O 



I^e2| 2 


Ks\ 2 


\Ves\ 2 


\V„i\ 2 + 






0.3 


0.5 


0.01 


0.11 
0.20 


0.94 
1.04 


0.95 
1.06 


0.3 


0.5 


0.005 


0.12 
0.18 


0.96 
1.03 
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fill 



Different values of R ejl and R eT for different combinations of 4 independant \ Vi 3 1 
and | V^x|+ span the range allowed by unitarity and positivity of | J| 2 ; these extremum 
values of (V^il 2 correspond to vanishingly small CP- violation. 
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Values of cos(S) for different combinations of parameters: cos(S) c is obtained from the 
central values of the parameters while cos(5)- and cos(5)+ are the extremum values 
allowed by the uncertainties (shown in percent, in parentheses next to each parameter). 
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Values of cos(8 e ff.) for different combinations of parameters in a non-radiative decay sce- 
nario: cos(5 e ff.) c is obtained from the central values of the parameters while cos(5 e //.)- 
and cos(8 e ff.)+ are the extremum values allowed by the uncertainties (shown in percent, 
in parentheses next to each parameter). 
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